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Abstract 

Partially depleted cores, as measured by core-Sersic model “break radii”, are typically tens to a few 
hundred parsecs in size. Here we investigate the unusually large (i? 7 /=o .5 = 4.57 kpc) depleted core 
recently reported for Holm 15A, the brightest cluster galaxy of Abell 85. We model the ID light 
profile, and also the 2D image (using Galfit-Corsair, a tool for fitting the core-Sersic model in 
2D). We find good agreement between the ID and 2D analyses, with minor discrepancies attributable 
to intrinsic ellipticity gradients. We show that a simple Sersic profile (with a low index n and no 
depleted core) plus the known outer exponential “halo” provide a good description of the stellar 
distribution. We caution that while almost every galaxy light profile will have a radius where the 
negative logarithmic slope of the intensity profile 7 / equals 0.5, this alone does not imply the presence 
of a partially depleted core within this radius. 

Subject headings: keyword: galaxies: elliptical and lenticular, cD — galaxies: individual (Holm I5A) 
— galaxies: photometry — galaxies: structure 


1 . INTRODUGTION 


Luminous early-type galaxies {Mb ^ —20.5 ± 0.75 mag) 
typic ally possess a core that i s partially depleted of stars 
(e.g. iDullo fc Grahanj 12014 and references therein). 
This is evident by a marked flattening of th e ir in¬ 
ner light distribution ([ Cing_&^inkowsO 119721 : 

IKormend^ 119821 : ILaueiHl983t IBvun et al.l Il996fl . The 
light profiles of these bright galaxies typically have a 


negative, inner logarithm ic slope 5,0.3, an d have been 
ht usin g the King mod e l (lKiiidll962L [196611. the “ Nuker 
law” (iGrillmair et al.l 11994 : _ Kormfind^eta]J_11994 


Lauer et al.lll995t iRavindranath et al.l 20011 Rest e t al.l 
2 OOIII. and the cor e -Sersic model (Graham et al.l 1200,4 
Truiillo et al.l 12004 iFerrarese et al. [200611 . The Nuker 
model was designed to investigate the inner-most regions 
of nearby early-type galaxies, and its double power-law 
nature was never intended to adapt to the full, outer, 


intrinsically curved i?^/”-like profiles of galaxies (e.g. 
iFaber et al.lIlQQTl their Section 2). The outer cur¬ 
vature in galaxies results in a profile whose slope changes 
as a function of radius. As a consequence, the slope of a 
htted Nuker model’s outer power law (/3) varies with the 
htted radial extent of a galaxy. Due to parameter cou¬ 
pling, it follows that all of the Nuker mo del parameters 
are a function of the fitted radial extent (|Graham et al.l 
120031 Figures 2-4). The coupling is such that the Nuker 
model “break radii” are heavily over-estimated, relative 
to the radius of maximum curvature in the actual pro- 
hle, i.e. the “break”, an d increasingly so a s the fitted 
radial exte nt is increased (iTruiillo et al.ll20'0^ . This con¬ 
tributed to lLauer et al.l (|2007[1 adopting the radius where 
the Nuker model has a negativ e logarithmic slope of 7 / 
equal to 0.5 (|Carollo et al.lll99^ as a measure of the core 
size. All galaxy light profiles, even those with no depleted 
cores , have a radius where 7 / equals 0.5. 

In iGraham et al.l (l2003f l it was shown how a Nuker 
model (with an inner power-law slope <0.3) can ap¬ 


proximate a Sersic profile (without any depleted core) 
if the profile has a lo w Sersic index and thus a shal¬ 
low inner profile slope. IDullo fc Grahaml (1201311 showed 
that this has occurred when mod elling several real galax - 
ies. Most recently, for example, iKrainovic et al.l (|2013ll 
report that their ATLAS^^ galaxy NGG 4473 is a 
core galaxy according to their Nuker model fit. How¬ 
ever, it actually contains an add itional nuclear compo - 
nent r ather than a deple t ed co re (|Krainovic et al.ll2006f) . 
which [Pullo fc Grahaml (120141 have shown through the 
use of the Sersic and core-Sersic models. The core-Sersic 
model was introduced, in part, to prevent confusion when 
a profile has a shallow inner slope but no central deficit 
of stars, and to provide more robust, physically meaning¬ 
ful radii, slopes, and flux deficits for partially depleted 
cores. 


Core sizes, as measured by the break radii of the core- 
Sersic model (Rh), are typically t ens to a few hundred 


parsecs (e.e. iTruiil 

o et al.l 12004 IFerrarese et al.l 

2006; 

Richings et al.ll2011 

: IDullo V Graba,mll2012ll2013l 

2014; 

Rush et al.l 1201311. 

These core-Sersic break radii are 


where the core-Sersic model has its ma ximum curvature, 
and, f or galaxies with depleted cores, IDullo fc Grahaml 
(j2012L their Section 5.2) revealed that this also matches 
well with the cusp radius^ Rj, where the negative log¬ 
arithmic slope of the projected light profile equals 0.5. 
As noted, galaxies without partially depleted cores can 
also possess such a radius, and therefore one still needs 
to establish if there is an inner deficit relative to the 
outer profile. Indeed, many galaxies have an inner 
profile slope shallower than 0.5, and 0.3, but do not 
have a partially deplete d core (|Graham fc Guznidi]l2003l : 
IDullo fc Grahamll2013l their appendix A.2). A thorough 
review of galaxy light profiles can be found in iGrahanil 

^ The cusp radius R^/ is defined as the radius at which the 
negative logarit hmic slope of the in tensity profile 7 / equals a pre¬ 
specified value l|Carollo et al.lfl997ll . Hereafter we will use to 
actually indicate i?^/=o. 5 - 
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(IMl . 

One of the most accredited scenarios for core for¬ 
mation attributes the depletion to the scouring ac¬ 
tion of black hole (BH) binaries formed during dry 
galaxy merger events. The binary depletes the cen¬ 
tre of a galaxy by ejecting stars via three-body in¬ 
teractions which result in the o r bital decay of the bi¬ 
nary (e.g. iBeeelman et al.l 119801: lEbisuzaki et al.l 119911 : 
iMilosavlieyic fc Merrittll2001l: lMerritdl200'^ 

lMerrittr(|2006ll suggested that the stellar mass depleted 
via the binary BH scouring mechanism scales as 0.5NM,, 
where N is the effective number of major dry merg¬ 
ers which the galaxy experienced. However, the ob- 
servational result that Mdej is typic a lly 0.5-4 M, (e.g . 
GrahamI 12004 iFerrarese et al.l 120061: iHyde et al.l 1 200 A 


Rush et al.l 120131 iDullo fc Grahaml 1201411 . implies that 
the core galaxies characterized by the larger Mdej /M, 
ratio should have experienced up to 8 major dry mergers. 
Such a number of major mergers near or above this figure 
is excessive if compared to the merger rates derived from 
the observation of clos e massive galaxy pairs in the local 
z < 0 .8) Universe fe . g.lBeh et al.l[2?i0^iDe Prooris et al.l 


[z < U .g) Universe (e . g.Itjell et aLII^UPtitlUe rrooris et al.l 
2QQ7I; [Casteels_eGS] [2 0141). A second mistery is that 


Savorgnan fc Grahaml (j201^ have revealed that the large 


bulges with over-massive BHs (or low velocity disper¬ 
sions) in the M,~a diagram do not have larger M^ef /M, 
ratios as expected f rom d ry mergers and as suggested by 
IVolonteri fc Giottl (j2013ll . 

The depleted core radius is indicative of the cen¬ 
tral mass depletion experienced by the galaxy, mod¬ 
ulo the pre-existing central stellar density profile. Scal¬ 
ing relations exist between the final merged mass 
(M,) of the central super-massive blac k hole (SMBH) 
and both the ejecte d stellar mass (e.g. iGrahamI 12004 
Ferrarese et al.ll2000fl and the c ore radii fe.g.lLauer et al.l 


20071 IDuHo fc Grahaml [2fn% IRusli et af .1 I20m 


These studies have described the M,-(core radius) cor¬ 
relation using a log-linear relation, albeit with large un¬ 
certainties due to significant scatter in the data. The 
ongoing effort to characterize the high mass end of the 
relations is necessary to investigate whether a single slope 
is appropriate to describe the diagram, as is cur¬ 

rently assumed in the range 10® < M,/Mq < 10^°. A 
bend in the relation may, for example, indicate 

different regimes of e fficiency for the BH scou ring mech¬ 
anism. For instance, iKulkarni fc LoebI (|2012t) have sug¬ 
gested that multiple (i.e., more than 2) SMBH systems 
can significantly increase the effectiveness of the scouring 
activity, and generate cores with a mass deficit up to five 
times M,. 

A non-linear relation may also indicate that additional 
mechanisms other than binary BH scouring are operat¬ 
ing. For example, a “kicked” SMBH can create an en¬ 
hanced de pleted core by crossing the nucleus multiple 
times (e.g. iRedmount fc Reedll989t iMerritt et al.ll2004 

iBovlan-Kolchin et al.ll2004 iGualandris &: Merrittll2008f) . 

In this scenario, the SMBH is placed on a radial orbit, in¬ 
tersecting with the nucleus, after the recoil acquired upon 
its creation from the coalescence of a specially-oriented 
BH binary. Such a recoil is generated in response to the 
linear momentum carried away by the anisotropic emis¬ 
sion of gravitational waves. 

Another possibility f or creating large cores is the 
“stalled binary” model (|Goerdt et al.l[20lO) . which pro¬ 


poses that the scouring activity is performed by a cap¬ 
tured in-falling object. In this picture, a “perturber” spi¬ 
rals towards a galaxian center due to dynami c al fric tion. 
The “stalled binary” scenario of iGoerdt et al.l (j20ir)fl pre¬ 
dicts core radii up to ^3 kpc and Mdef scaling as the 
mass of the perturber. In doing so, it exerts a tidal ac¬ 
tion on the central mass distribution, shredding it and 
creating a partially depleted core. This model predicts 
that the core radius will correspond to the orbit at which 
the infall of the perturber stalls due to a reduced efh- 
ciency of the dynami cal friction. In yet another scenario, 
iMartizzi et al.l (|20I2t) argued that the feedback action of 
an active galactic nucleus (AGN) might be an impor¬ 
tant facto r in the creation of c ores. The AGN feedback 
models bv iMartizzi et ahl (|20I2f) easily produce core radii 
of sizes up to ~I0 kpc although not yet confirmed. In 
this model, AGN-driven gas outflows generate fluctua¬ 
tions in the gravitational potential of the central region, 
from which stars are removed during the subsequent re- 
virialization and adjustment process. The expulsion of 
gas from the inner region might also induce an adiabatic 
expansion of the central stellar distribution, hence flat¬ 
tening the central mass density profile. 

Based on the extrapolation of the relation between 
the luminosity of the host spheroid (L ) and M, (e.g. 
iGraham fc Scottl I20I3I : IGrahamI l20I5bll . some Bright¬ 
est Cluster Galaxies (BCGs) are expected to host ex¬ 
tremely massive BHs (M, > 10^® Mq). Recent N-body 
simulations including the effects of dark matter, black 
holes, as well as baryons, have shown that black hole 
scouring in BCGs can create cores as large as 3 kpc 
(jLaporte fc Whit ^12014( 1. Therefore, the study of BCGs 
is important to explore the scarcely-populated high-mass 
end of the M,-Rd diagram, and hence provide better con¬ 
straints on the Mt~Rb scaling relation and the formation 
physics involved. 

1.1. The case of Holm 15A 

iLopez-Gruz et all (j20l4 hereafter: LCI4) studied the 
BCG Holm ISA (D = 253 Mpc)^, located within the 
galaxy cluster Abell 85, and reported the discovery of 
the largest depleted core known using the Nuker model. 
They found a cusp radius R^ = 4.57 ± 0.06 kpc, which 
supersedes the record previo usly reported for the BCG 
in Abell 2261 {R^ ~ 3.2 kpc: iPostman et ani20I2h . The 
cusp radius of Holm ISA is not only mor e than three 
times larger than the average R^ obtained bv iLauer et al.1 
(|2007tl using the Nuker model for their sub-sample of ^60 
BCGs, but is also significantly larger than the biggest 
recorded co re-Sersic model bre ak radius (i?& 1.5 kpc) 

reported bv iHvde et al.1 l)2008ll for the massive elliptical 
galaxy SDSS J09I944.2-b56220I.I. 

LC14 derived R^ fitting a 2D Nuker model to an 
image of Holm I5A within a major-axis radius of 
^80 kpc. They confirmed their measurement of R^ non- 
parametrically from the ID radial light profile (i.e., ap¬ 
plying the definition given in Footnote [T]). However, 
their inner light deficit (Tde/) was not calculated using 
the Nuker model, but rather as the difference between a 

^ Luminosity distance from NED, relative to the Cosmic Mi¬ 
crowave Background, assuming Hq = 67.30 km s~^ Mpc~^, = 
0.315, Da = 0.685 (Plancfc-I-WMAP; IPlanck Collaboration e t al.l 
l20ll) . 
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double-Sersic fit (intended to represent the actual galaxy 
light distribution) and a ide VaucoulenrsI (jl948[l pro¬ 
file (intended to reproduce the “pristine” light profile 
prior to the redistribution of the inner core light). This 
unconventional light deficit was the n used to derive a BH 
mass using the f-M, relation of iKormendv fc Bended 
(l200l . giving an exceptionally high M, ~ 10^^ Mq. 
LC14 assessed the viability of this value by comparing it 
against the value of M, estimated using other methods. 
Since no direct (dynamical) measurement of the mass of 
the SMBH of Holm 15A is available, LC14 resorted to 
using the scaling relations between M, and: the stellar 
velocity dispersion (cr); the total lumin osity of the bulge; 
the N uker model break radius; and (|Kormendv &: Hoi 
I2013L and references therein). This analysis ultimately 
lead them to favor a more conservative M, ^ 10^° Mq. 

Is the exceptionally large core of Holm 15A really due 
to a deficit of light relative to the inward extrapolation of 
its outer light profile? Is there an obvious and dramatic 
downward bend to its inner light profile? We have further 
investigated the case of Holm 15A by performing a de¬ 
tailed analysis of the light distribution of the galaxy. For 
the first time, we apply the core-Sersic model and check 
if a depleted core is warranted over a core-less model (i.e. 
a Sersic model). 

This paper is structured as follows. In !j2l we present 
the data and the procedure used for modeling the ID 
light profile and the 2D image of the galaxy. In (j3] we 
show that the light in Holm 15A is well fit by a three- 
parameter Sersic function plus an exponential “halo” and 
as such it does not appear to have a partially depleted 
core with a well-defined break in the light profile. We 
summarize our conclusions in (J4l 

2 . DATA 

We have used an r-band^ image from the wide-held 
MegaPrime camera mounted on the Canadian-French- 
Hawaiian-Telescope (CFHT). The image was retrieved 
from the Canadian Astronomy Data Center'^ (CADC). 
This enabled us to make a direct comparison with LC14, 
who used Kitt Peak National Observatory (KPNO) and 
CFHT data, in the R and r bands, respectively. We 
have based our analysis on the CFHT data because — 
to our knowledge — they represent the publicly avail¬ 
able images for Holm 15A with the best seeing quality. 
A narrow point spread function (PSF) is important to 
detect core structures, whose projected sizes, even in the 
closest galaxies, are of the order of a few arcseconds or 
less. 

The detector of the CFHT MegaPrime camera is com¬ 
posed of 36 CCDs (2112x4644 pixels each) covering a 
sky area of ^1 square degree. Every MegaPrime im- 
age in the CADC is reduc ed via the Elixir pipeline 
(|Magnier fc Cuillandrel 1^0041) . which performs the basic 
data calibration, including bias subtraction, flat-helding, 
de-fringing, and astrometric calibration. The CFHT 
data used by LC14 was a single 120 s exposure. We 
searched the CADC archive for a similar data set, and 

® Nominally, in the r.MP9601 filter. This band corresponds to 
the SDSS r-band, to which we will refer hereafter, 
http://www2.cadc-ccda.hia-iha.nrc-cnrc.gc.ca 


discovered that all the images were severely affected by 
scattered light. This scattered light, which appears as 
a diffuse radial gradient across the image, is known to 
plague the MegaPrime data and it is due to unwanted 
internal reflections through the optics. This wavelength- 
dependent issue, especi ally promin e nt in the r-band, was 
extensively reviewed by iDuc et all (1201511 . 

Due to the lack of an analytical description for the 
scattered light, it could not be readily disentangled from 
the smoother sky background. We therefore modeled 
the combined background (sky -|- scattered light), and 
subtracted the contribution from the two components 
at the same time. Since Holm 15A easily fits within 
one single MegaPrime CCD, still allowing ample mar¬ 
gins to determine the background, we performed this op¬ 
eration and the other analyses only on the CCD which 
was centered on the galaxy. First, we ran S Extrac¬ 
tor (iBertin fc Arnouta 1 19961 1 to determine the i/go (i.e. 
the radius encompassing 90% of the flux) of each source 
in the image, and to obtain a median estimate of the 
remaining background flux and of the background root- 
mean-square (RMS) fluctuation. Each detected source 
was masked using an elliptical region with a major axis 
equal to twice the size of the i/go of the source, and with 
axis ratio and position angle (P-A.) as determined by 
SExtractor; in addition, we masked any ba fluctua¬ 
tion abo ut the median back ground flux. Finally, we used 
Galfit (iPeng et al.l 12010( 1 to ht a 2D gradient to the 
background, adopting the SExtractor median back¬ 
ground as the initial guess for the yalue at the image cen¬ 
ter. Eigure [T] shows our background-subtracted r-band 
image (see Table [T] for image specifications). 

2.1. The ID radial light profile 

An object mask was created using the SExtractor de¬ 
tections, and then further refined by hand after yisual 
inspection to exclude cosmic rays, hot pixels, and smaller 
sources oyer-lapping with Holm 15A. The right panel of 
Figure [1] shows the CFHT mosaic with the masked areas 
down-scaled by an arbitrary amount; we also oyerplot the 
elliptical isophotes deriyed using the IRAF. ellipse task 
(|Jedrzeiewskilll98^ . 

Figure [2] shows the radial profiles along the semi-major 
axis (SMA) obtained from IKAF.ellipse for the: r-band 
surface brightness {fir)', ellipticity (e); Ath harmonic de- 
yiation from perfect ellipticity {BA ’’boxiness/diskiness” 
parameter); P.A.; and isophote centroid shift. We ob¬ 
serve that beyond ^^2", the ellipticity is steadily increas¬ 
ing, except for a flattening in the range 10"-30". Increas¬ 
ing ellipticities have been shown to be common for BCGs, 
and they can be due to the projection of the prolate (or 
triaxial) structure of the outer regions (typically b eyond 
^45 kpc) of these objects (e.g. iPorter et al.lll991r). 

The fitting of models to the light profile (see ii2.3l) haye 
been performed yia a Leyenb erg-Marquardt minimiza - 
tion procedure, as described in iDullo fc GrahamI (|2014[) . 
We chose to limit the fit to fir < 25.5 mag arcsec“^ 
(approximately corresponding to SMA < 90 kpc), based 
on where the residual large scale background gradient 
(i.e., after the background subtraction process described 
aboye) becomes comparable to the galaxy light gradient. 
This limit roughly corresponds to the radius at which the 
elliptical isophotes stop being concentric (see Figure [U 
right panel). 
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Fig. 1.— CFHT-MegaPrime r-band image (left) and relevant mask (right) for Holm 15A. 

The masked areas (see have been arbitrarily decreased in intensity so as to still show the contaminating objects. We overplot the 
elliptical isophotes from IRAF. ellipse (the boxiness of the ellipses is not considered in this representation). The dashed ellipse corresponds 
to the limiting surface brightness at which we truncated our ID analysis {fir 25.5 mag arcsec“^; see This ellipse also corresponds 
to the physical extent of the 2D fit (everything outside the dashed curve was masked for the 2D fit; see [[2}. 


TABLE 1 

CFHT-MegaPipe image characteristics 


Target 

(1) 

RA (J2000) 
[hh:mm:ss] 

(2) 

Dec (J2000) 
[dd:mm:ss] 

(3) 

D 

[Mpc] 

(4) 

m — M 
[mag] 

(5) 

Camera/Filter 

(e) 

Exposure 

[sec] 

(7) 

Scale 

[^Ypixel] 

(8) 

Holm 15A 

Oo'i41"»5os 5 

-09°18'11" 

253 

37.02 

MegaPrime/r 

120 

0.186 


Note. — Details of the CFHT-MegaPrime image used for the current work. 

Target name. Target coordinates from NED. Luminosity distance from NED, 

corresponding to a redshift z ~ 0.057 (see Footnote il). Distance modulus. CFHT 
camera and filter. Total exposure time. Image pixel scale. 


2.2. The 2D image 

Modeling of the 2D ligh t distribution was performed us¬ 
ing Galfit-Corsair® ()Bonfinil[2014D . a tool developed 
to include the core-Sersic model into the Galfit fitting 
algorithm®. 

To maximize the consistency with the ID analysis, we 
restricted the fit within the Hr ^ 25.5 mag arcsec“^ 

® Iwww.astronomy.swin.edu.au/'pbonfini/galfit-corsair/l 

® Galfit-Corsair is retro-compatible with Galfit. 


isophote (see H2.1|l by masking every pixel ontside it (see 
Figure m right panel). The Galfit-Corsair mask was 
hence obtained by combining this filter with the object 
mask des cribed before . The PSF template was built with 
PSFEx (|Bertinll2011l) . combining ^20 stars with signal- 
to-noise ratios {S/N) > 100. Since the variation of the 
PSF FWHM was negligible across the CCD (less than 
5%), we selected stars from all areas. The FWHM of the 
resulting PSF model is ~0".75, in agreement with the 
value measured for real stars in the image. Finally, the 
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Fig. 2.— IHAF.ellipse major-axis radial profiles for: r-band surface brightness (top-left), ellipticity (top-right), Ath harmonic deviation 
from perfect ellipticity (BA; middle-left), position angle (middle-right), and isophote centroid shift along the x-axis (bottom-left) and y-axis 
(bottom-right) with respect to the innermost isophote. 


weight (“sigma”) image was constructed using the inter¬ 
nal Galfit algorithm, to which we supplied an estimate 
of the sky RMS as measured on the image by using SEx- 
TRACTOR. 

We warn the reader that the varying ellipticity of 
Holm 15A (Figure [2l top-right panel) can pose a chal¬ 
lenge for the 2D fitting, because in Galfit-Corsair, as 
with Galfit, each galaxy component is associated with 
a single ellipticity and position angle. When the compo¬ 
nents of real galaxies have ellipticity gradients, the model 
does not allow for this. That is, even when a galaxy is de¬ 
composed into multiple components, the problem might 
persist whenever the ellipticity profile has a significant 
gradient over the range where any single component dom¬ 
inates. This issue is of no concern when modeling the ID 
profile, which is simply extracted along the major axis. 
As such, small discrepancies are expected between the 
ID and 2D analyses. 


2.3. The choice of models 

To investigate the presence of a partially depleted core 
in Holm 15A, we separately fit a seeing-convolved Sersic 
model and a seeing-convolved core-Sersic model, in both 
ID and 2D. As noted and modelled by LC14, Holm 15A 
has a halo of light around it. The “bump” in the radial 
light profile at around 35 kpc, which also corresponds to 
a step in the ellipticity profile (see Figure[2|), does indeed 
suggest that the galaxy hosts a second component. We 
therefore added an outer [seeing-convolved] exponential 
function to the Sersic/core-Sersic “bulge”, capturing the 
“halo” of Intra-Cluster Lig ht (ICL). An e xpone ntial halo 
model was used following iSeigar et al.l (I2007D . who re¬ 
vealed t hat the halos of BC Gs are typically exponential; 
see also lPierini et al.l (|2008f) . We did additionally model 
the ICL with a Sersic model, but we found that it 
yielded a Sersic index of ^1, i.e. an exponential model. 

The best-fit parameters obtained for the aforemen- 
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tioned models are reported in Table [H while the model 
profiles are represented in Figure[3](lD fit) and in the left 
and central panels of Figured (2D fit). We stress that 
the profiles shown in Figure [4] represent a ID projection 
of the 2D images and 2D models rather than what was 
actually minimized, namely the difference between the 
2D image and model, as shown in lower panels of Figure 
[4l This projection was performed measuring the surface 
brightness of these images along the isopho tes i dentified 
with \JiAY.ellipse in our ID analysis (see 42.11) . In the 
next section, we will first mention the small discrepancies 
between the results of the ID and 2D modeling, and then 
proceed to discuss the selection of the best-fit model. 

2.4. Fit comparison and choice of best-fit model 

Overall, our results from the modeling of the ID light 
profile and the 2D image agree very well (see Table El). 
However, the ID projection of the 2D residuals (Figurell]) 
show a more pronounced pattern than seen the ID resid¬ 
uals (FigureEl). This is partly an artifact of the particular 
choice of projection, and partly a direct consequence of 
each 2D model component of Galfit-Corsair having a 
single center and ellipticity weighted over the whole ex¬ 
tent of the fit area. Here we will review the implications 
of these caveats over two specific regions. 

The innermost '^0" .5. — In the projection shown in Fig¬ 
ure IH the 2D model seems to significantly underestimate 
the light at radii smaller than '^0".5. Within ^0".5, 
the isophote centers identified by IRAF. ellipse — which 
can trace the centroid of the real galaxy at the sub-pixel 
level — drastically shift (see Figure El bottom panels). 
When the photometry is performed over the 2D model 
along the same isophotes, the centroid of the Galfit- 
CORSAIR model are offset, hence the model brightness is 
underestimated within the ~0".5 region. 

However, since a corresponding tiny “bump” is visible 
in the ID surface brightness profile (Figure El top left 
panel) at those radii, we investigated whether this small 
feature should be attributed to an intrinsically lopsided 
light distribution, or instead to the presence of a mis¬ 
aligned point-source. If there was in fact a point-source, 
it could be an unresolved nuclear star cluster or an offset 
AGN. Indeed, LGI4 found evidence for AGN emission in 
the optical spectrum of Holm 15A, which might be visi¬ 
ble in our r-band image. We tried to model this potential 
inner feature by adding a point-source to the 2D mod¬ 
els (nominally, a PSF function), but we could not obtain 
a successful fit due to the low surface brightness of this 
potential component relative to the surrounding galaxy. 
We also tried to perform a similar fit in ID using a PSF 
and a PSF-convolved Gaussian (i.e., a slightly extended 
source), although our ID residuals also did not warrant 
such an additional component. From this fit we obtained 
a PSF-convolved Gaussian with FWHM comparable to 
the PSF FWHM (as expected for a point-source), and 
an integrated magnitude for any potential nuclear com¬ 
ponent of ^24.8 mag, which we consider as an upper 
limit. 

The 1"-j" region. — As discussed above, the elliptic¬ 
ity profile of Holm ISA (Figure El top-right) suggests 
the presence of two galaxy components (i.e., spheroid 
and halo). Assuming that the 10"-30" interval marks 
the transition between them, our 2D fit seems to recover 


the ellipticity at the outer edge of each component (see 
Table El)i most probably due to the higher number of 
pixels along larger isophotes. The single ellipticity of 
the Galfit-Corsair components is the cause for the 
pattern visible in the center of the residual image (Fig¬ 
ure IH bottom-left), which appears as a I"-4" structure 
oriented in the direction perpendicular to the galaxian 
major-axis (i.e., along the NE-SW direction). This hap¬ 
pens because in that region the image isophotes have a 
much more circular shape ((e) O.I) than the 2D Sersic 
component (for which we obtained e = 0.2). In the ID 
projection of the 2D residuals (Figure IH top-left), this 
pattern manifests as a sinusoidal feature. 

As a consequence, the ID and 2D core-Sersic-l— 
exponential fits show some differences in their inner re¬ 
gions. The ID model shown in right-hand panels of Fig- 
ureElyields a small (sub-pixel) break radius {Rb ^ 0".14), 
generally consistent with being a core-less galaxy^. Ev¬ 
ery core-Sersic fit can effectively reproduce an intrinsic 
Sersic profile by sufficiently minimizing the break radius 
of the core. Although a core is partially resolved in the 
2D model (middle panels of Figured]), it does not seem 
to improve the fit significantly. In particular, we observe 
that the I"-4" 2D residual feature described above is still 
present, and — although slightly reduced — still shows 
the same pattern, hence reducing the likelihood that it 
is related to a depleted a core. 

On the other hand, outside the problematic in¬ 
ner ^0".5 region, the Sersic-|-exponential and core- 
Sersic-{-exponential models fit the data comparably well, 
yielding similar residuals (see Table El and Figures |3| and 
lH). What has occurred is that the core-Sersic core param¬ 
eters were driven by the fit algorithms to partially com¬ 
pensate for the effects of the ~0".5 “excess” (ID case) or 
the varying ellipticity (2D case), rather than accounting 
for an actual depleted core. 

Persuaded that the simpler Sersic-|-exponential model 
provides the most appropriate description, we decided 
to refine the 2D Sersic-|-exponential fit using nested 
Galfit-Corsair models to see if we could approximate 
the varying ellipticity of Holm 15 A. In addition to an el¬ 
lipticity gradient, the radial range corresponding to the 
inner Sersic component is also associated with an irreg¬ 
ular position angle profile which is wildly varying within 
the innermost 4" (while it is remarkably constant in the 
range where the exponential component dominates; see 
Figure El middle-right). Therefore, we paid special at¬ 
tention to improve the fit of the inner component, and 
we constructed our refined model by using two Sersic 
components with different axis ratios to describe the 
“bulge”, plus the previously used exponential function 
for the “halo” (FigurejH right panels). This extra central 
component represents a corrective factor, rather than a 
distinct element of the galaxy. The parameters obtained 
with this new fit are reported in Table El The best-fit ex¬ 
ponential halo component is practically unchanged with 
the addition of this new inner Sersic component, indi¬ 
cating that the correction indeed acted mostly on the 
“bulge”. Moreover, the corrective component has a lu¬ 
minosity one order of magnitude fainter than the other 
components, so that the predominant Sersic component 

^ Recall that the image has a seeing of ~0”.75. 
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Fig. 3.— ID light profile analysis: results of the Sersic+exponential (left), and core-Sersic+exponential (right) model fit to the major- 
axis, r-band surface brightness profile of Holm 15A. 

In the left panel, the red curve shows the Sersic component. In the right panel, it represents the Sersic portion of the core-Sersic component. 
The exponential function which dominates at large radii is indicated by the violet short-dashed curves. The solid curve represents the 
complete fit to the profiles, with the RMS residuals, A, about each fit given in the lower panels. Note that the core part of the core-Sersic 
model fits an apparent slight “excess” of light over the inner 0".5, rather than a light “deficit”. 


still resembles the previously obtained “bulge” (including 
its position angle and aspect ratio). With this refinement 
the 2D fit drastically improved, now yielding fiat resid¬ 
uals over all the fitted range beyond ~0".5 (see Figure 
m bottom-right panel). The apparent (false) excess seen 
here, in the ID profile is a result of the centroid shift 
issue discussed earlier. 

3. DISCUSSION 

Extended halos around BCGs are commonly observed. 
Their presence is not only indicated by the need for 
adding an “envelope” component when fitting models, 
but also by the twist in the outer isophotes of BCGs, 
which has been interp reted as a signature of galaxy accre¬ 
tion/ interaction (e.g. iPorter et aLlIlQQll : iGonzalez et al.l 
l2005f) . This picture is supported by numerical simula¬ 
tions, which showed that BCGs, being at the centers of 
large dark matter haloes, underwent a vigorous history of 
both minor and major merging (e.g. iDe Lucia fc Blaizotl 
[2003), and hence may be surrounded b y a halo of tidally 
stripped stars te.g. lGooner et al.]l2014[ l. 

In fa ct, in their survey of nearby BGGs. lDonzelli et al.l 
dmi herafter DMIl) had already reported that 
a Sersic-bexponential model fits the light profile of 
Holm ISA. We improve on this by: (I) investigating 
the presence of a depleted core (DMII did not per¬ 
form this test plus their data had a FWHM 1"- 
2" and their fits were limited to a minimum radius of 
1.5XFWHM); and (2) reaching deeper in limiting sur¬ 
face brightness ('^25.5 mag arcsec”^ compared to their 
24.5 mag arcsec”^). DMII thus limited their analysis to 


^60", while our fit extends up to ^80", hence providing 
a better constraint on the extended galaxy halo light. 
Our ID Sersic-bexponential fit agrees well with the simi¬ 
lar fit by D MII. Given the luminosity~n relation and the 
findings in iGrahamI (II996I1 , it may be surprising to note 
that the Sersic index that we derive for the main “bulge” 
of a bright galaxy like Holm I5A {My ^ -23.8 mag) is 
so small (n 1). However, this low Sersic index can 
be compared to the low Sersic indexes measured for the 
inner components of several BCGs (e.g. DMII), as well 
as in a couple of c D galaxies (NGC 4874 and UGC 9799; 
iSeigar et al.ll2007ll . These low Sersic indexes might be a 
result of a dramatic galaxy re-shaping due to one/some of 
the processes mentioned in the Introduction. Although 
these low-n profiles have rather flat inner regions, we are 
not dealing with a profile that displays a clear break and 
downward deviation at small radii — as observed with 
the traditional partially depleted cores in galaxies with 
larger (n > 3) Sersic indexes. 

In Figure El we display the in n er ga laxy components 
from the fits by iDonzelli et al.l (1201111 (their Table 2; 
dashed lines) and ([Seigar et al.! 120071 their Table 2; 
magenta solid lines), after selecting objects similar to 
Holm 15A, i.e. with Sersic-l-exponential decomposition, 
and hayi ng inner sphero i ds wi th Sersic n < 1.5. We high¬ 
light the IDonzelli et afl (j2011t) model for Holm 15 A with 
a black dashed line, which can be obseryed to closely 
match the profile we derived (solid green line). All 
the surface brightness profiles are expressed relative to 
the R-band, and, when necessary, have been converted 
from the observed bands using the colour conversions 
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Fig. 4.— 2D image analysis: results of the fits to the 2D image of Holm 15A performed assuming: a Sersic+exponential model (left); a 
core-Sersic+exponential model (middle); and a Sersic+exponential model plus an additional inner Sersic component intended to compensate 
for the single ellipticity of the 2D model components (right). 

Top.- The green data points represent the major-axis surface brightness profile measured over the isophotes defined using IKAF.ellipse 
(i.e. the same measurement presented in Figure [JJ. The curves represent the surface brightness profile of the model images measured over 
exactly the same isophotes. The continuous curves show the global models, while the dashed curves represent their sub-components. We 
stress that these are not fits to the ID profile, but rather surface brightness measurements (projections) of the 2D models. The pentagon 
indicates the location of the core-Sersic model’s break radius. The panels underneath the profiles represent the data residuals about the 
fitted models, first expressed in terms of the difference in surface brightness, and then in terms of residuals (in units of counts) divided by 
the standard deviation as measured on the “sigma” image. Bottom.- The actual residual images that were minimized by Galfit-CorSAIR. 
Masked objects are highlighted as in Figure^ 














































TABLE 2 
Fit Results 


Model / Component 

Mr 

'^comp,r 

Rb 

Rb 

a 

7 

Re 

Re 

n 

e 

RA. 


+0 

<0 



[mag/arcsec^] 

[mag] 

[kpc] 

[arcsec] 



[kpc] 

[arcsec] 



[deg] 

[mag/arcsec^j 

[mag] 

[mag] 

(1) 


(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

(12) 

(13) 

(14) 

(15) 

ID 

Sersic+exp. | 

Sersic 

exponential 

21.60 

21.93 

13.78 

13.32 





13.09 

46.49 

11.87 

42.15 

1.0 

[1.0] 



} 0.015 

12.50 

-24.52 

core-Sersic+exp. | 

core-Sersic 

exponential 

19.83 

21.94 

13.79 

13.33 

0.15 

0.14 

[2.0] 

0.28 

12.86 

46.60 

11.66 

42.25 

0.9 

[1.0] 



} 0.014 

12.51 

-24.49 

2D 

Sersic+exp. | 

Sersic 

exponential 

21.59 

21.83 

13.82 

13.47 





12.81 

44.08 

11.61 

39.96 

1.0 

[1.0] 

0.19 

0.38 

143.9 

148.0 

} 0.119 

12.61 

-24.41 

core-Sersic+exp. | 

core-Sersic 

19.96 

13.76 

0.81 

0.73 

1.2 

-0.12 

13.30 

12.06 

1.0 

0.20 

144.0 

} 0.119 

12.62 

-24.40 

exponential 

21.99 

13.53 


46.90 

42.52 

[1.0] 

0.40 

148.2 

r 

Sersic 

21.62 

13.85 





13.19 

11.96 

0.9 

0.21 

143.9 

I 



Sersic+exp.+Sersic < 

exponential 

21.85 

13.46 





44.37 

40.23 

[1.0] 

0.37 

148.2 

> 0.118 

12.61 

-24.41 

Sersic (corr.) 

22.67 

17.82 





4.13 

3.74 

0.3 

0.24 

54.9 




Note. — Best-fit parameters from our ID and 2D analysis. Missing values are not relevant to the model/component under consideration. 

Fit model/component. r-band surface brightness at the: break radius (core-Sersic), effective radius (Sersic), or at i? = 0 (exponential). Total r-band 

magnitude of the component. For the ID models, we integrated the surface brightness profile assuming a constant ellipticity (e), which we chose as the ellipticity 
at the component’s effective radius ((e) ~ 0.2 for the inner components, (e) ~ 0.3 for the outer exponential components; see Figure [2)). Break radius in units of 
kilo-parsecs (at the distance of Holm 15A, l"= 1.103 kpc). Break radius in units of arcseconds. Alpha parameter for the core-Sersic model; this parameter 
has been held fixed (i.e. not solved for as a free parameter) in the ID analysis. Inner power-law index for the core-Sersic model. Effective radius (Re) of the 
component in units of kiloparsecs; the Re of the core-Sersic profile is defined to be the effective radius of the whole model in the ID analysis, while it represents the 
Re of the Sersic part of the composite model in Galfit-Corsair; since in this case the break radius turns out to be very small, the two effectively coincide. 
Effective radius in units of arcseconds. Sersic index (for n = 1 the Sersic model corresponds to the exponential model). Component ellipticity. Component 
position angle (North = 0°). Dispersion of fit residuals. For the ID fits we adopted the RMS of the difference between the light profile and the model values 
(see Figure [3|. For the 2D fits we adopted the inner 68% (±lcr) interval of the distribution of values from the residual image pixels (after applying the same masking 
used for the fit). Note that the former yields numerically smaller values than the latter because the residual data points are averaged over the whole isophote. 
Rest-frame, extinction-correct ed r-band magn itude of the model. A 51og(l-|-z) 0.12 mag magnitude dimming was applied using the distance reported in Table [T] 

Following the prescriptions of IPoggiantil (Il997f) for early-type galaxies, we also applied an r-band A-correction Kr = 0.06 mag. We used the galactic extinction from 
NED (~0.09 mag). Absolute, rest-frame, extinction-corrected r-band magnitude of the model (distance modulus is provided in Table [T]). 
t Values refer to the SDSS r-band filter, ABMAG system. The zero-point for the calibration is provided in the CADC image header. 


CO 
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Fig. 5.— Inner components of local BCGs with low Sersic in¬ 
dexes. The dashed lines (magenta solid lines) repre sent the Sersic 
compo n ents of the Ser^c H-exponential fits from IDonzelli et al.l 
J2011])_JSei^aT et_ajJ| [^07l) . The black dashed line represents the 
IDonzelli et al.l J) inner component of Holm 15A, while the solid 
green line shows its spheroidal Sersic component from our analysis. 
All the surface brightness profiles are expressed relative to the R- 
band, and, when necessary, have been converted from t he ob served 
bands using the colour conversions by IFukugita et al.l 11199511 . 


bv IFukugita et all l)1995h . Holm 15A does not seem to 
show any striking peculiarity in this representation, ex¬ 
cept that it has the l owest central surface brightness in 
the [sub-]sample from IDonzelli et al] (|2011h . This obser¬ 
vation strengthens the idea that Holm 15A is an other¬ 
wise somewhat standard representative of the population 
of local BCGs with low Sersic index spheroids. 

While LC14 favor the Nuker law for their descrip¬ 
tion of the central regions of Holm 15A, they also use a 
double-Sersic fit (comparable to our Sersic-fexponential 
model) to recover the total galaxy light, obtaining 
mv = 13.28 mag. We can compare this value with 
the total magnitude obtained from our 2D Sersic- 
-l-exponential model, before rest-frame and extinction 
correction (m^ = 12.88 mag; see Table [5] and table 
notes). We converted our magnitude, defined in the 
SDSS r-band (AB system), to the Johnson/Cousin V- 
band magnitude in the Vegamag system (my) by ap¬ 
plying a {V — r') = 0.36 ma g color correction, as t ypical 
for early-type galaxies le.g. IFukugita et HI 119951) . We 
obtained my = 13.24 mag, comparable to the measure¬ 
ment of LCi4. 

An overview on how we came to believe in massive 
black holes, and the development of the scaling relations 
associat ed with t h em, ca n be found in the extensive re¬ 
view bv iGrahamI (l2015ari . LG 14 estimated the mass of 
the BH in Holm 15A using several scaling relations (their 
Table 2), and concluded that the best estimate for M, 
is ^10^^Mq. Of the relations they used, the only one 
which directly connects to the assumed m orphological 
profile was the M,-Rb relation taken from iRusli et al.l 
(|2013ll . from which they obtained M, ~ 1.7x IO^^Mq us¬ 
ing their Nuker law break radius. This is an o rder of mag¬ 
nitude above their preferred estimate. The IRusli et HI 
(|2013f l scaling relation was however constructed using the 
core-Sersic break radius {Rb^), which is smaller than the 


Nuker break radius (R^) for the reasons discussed in the 
Introduction. However, we found no convincing evidence 
in Holm 15A for a depleted core relative to the inward 
extrapolation of the outer light profile and as such we do 
not deem it appropriate to use the relation. 

In the binary black hole scouring scenario, the core 
mas s deficit relate s to the mass of the final merged SMBH 
fe.g. lMerrittll2006ll . The “depleted” mass can be inferred 
from the luminosity deficit with respect to a model rep¬ 
resenting the “unperturbed” galaxy. In LG14, they use 
a de Vaucouleurs profile as their original, unperturbed 
model to compare against their double-Sersic fit. From 
a visual inspection, their de Vaucouleurs profile matches 
the data over the range of galactocentric radii: 10-40 kpc 
(LC14; their Figure 2). However, given that galaxies and 
BCGs exhibit a variety of Sersic indexes (e.g. DM11), it 
might not be that the original (i.e., pre-depletion) profile 
of the galaxy had n = 4 (corresponding to a de Vau¬ 
couleurs profile). Indeed, the application of a de Vau¬ 
couleurs model appears to have produced an artificial 
over-sized luminosity deficit. 

The discrepancies between our results and those of 
LC14 are related to the adopted paradigm for the de¬ 
scription of the surface brightness distribution, i.e. the 
core-Sersic or the Nuker framework. There is no formal 
mistake in the analysis of LC14, and in fact we were 
able to reproduce their Nuker parameters when fitting a 
Nuker model in 2D with Galfit. However, for the rea¬ 
sons highlighted in the Introduction — and especially for 
its ability to discern real cores from the flat inner slopes 
of low-n Sersic profiles — we favor the core-Sersic model, 
which found no evidence for a depleted core larger than 
the resolution limit. 

Finally, we remark that the depletion of stars from the 
core of a galaxy due to coalescing black holes, and other 
mechanisms, preferentially removes stars on radial or- 
bits, leaving an exce ss of stars on tangential orbits (e.g. 
iThomas et al.l[^014f) . Integral-field kinematic data may 
therefore be of benefit to help identify or reject the pres¬ 
ence of a core depleted of stars. 

4. CONCLUSIONS 

We performed a surface brightness analysis of the BCG 
Holm 15A using a CFHT-MegaPrime r-band image to 
investigate the claim that this galax y has the largest 
deple ted core ever detected (4.57 kpc: lL6Dez-Cruz et al.l 
[2OIII . 

We fit the ID light profile and the 2D image of 
Holm 15A to compare a core-less galaxy plus envelope 
model (Sersic-bexponential) against a core-Sersic galaxy 
plus envelope model (core-Sersic-1-exponential). We ob¬ 
tained good agreement among the best-fit parameters de¬ 
rived with the ID and 2D methods (see Table la, mod¬ 
ulo some minor differences predominantly attributable 
to ellipticity gradients f il2.4l) and a varying center with 
isophotal radius. In order to approximate the varying 
ellipticity of the inner galaxy, the 2D models — having 
components with fixed ellipticity — required the addition 
of a “corrective” component to the Sersic-l-exponential 
model. 

We find that the core-Sersic model does not provide 
an appropriate representation of the galaxy light distri¬ 
bution. In particular, in the ID description the core has 
a slight excess of light, while in the 2D model the inner 
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power-law of the core-Sersic fit does not represent an ac¬ 
tual real core, but rather compensates for the ellipticity 
gradient f 112.41) . We therefore conclude that the galaxy is 
core-less and we favor the idea that its light distribution 
is best described by a simple Sersic profile with a low in¬ 
dex n plus an exponential “halo” component as included 
by LC14. 
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